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Executive  Summary 


The  woik  rqwited  in  this  documoit  establishes  the  theoretical  feasibility  of  a  straightforward 
and  practical  spectroscopic  method  for  determining  the  temperature  of  ion  beam  target  plasmas 
derived  from  hydrocarbon  foams.  The  technique  is  most  qjprppriately  tq^lied  to  die  time  frame 
following  the  passage  of  the  ion  beam ,  when  direct  beam  excitation  of  the  atomic  states  of  the  target 
material  has  ceased.  For  die  temperature  range  of  50-200  eV,  it  is  found  diat  unifotmly  doping  die 
foam  widi  NaF  and/or  oxygen,  each  at  atomic  densities  of  about  6  x  10^'  cm~’,  will  result  in  die 
productitm  of  charactetisdc  K  shell  x-n^  features  from  diese  elements.  For  reasons  discussed  in 
dx  many  of  the  qiectral  lines  are  predicted  to  appear  in  absorpdrm.  The  temperatures  are 
primarily  differairi a»ed  by  the  presence  or  absence  of  various  lines  from  the  six  K  shell  ionizadoo 
stages  of  these  three  elements.  Some  of  die  lines  would  be  seen  in  emission  in  die  iqiper  part  of  die 
temperature  range  considered.  This  {dienmnenon  dso  is  useful  as  a  temperature  indicator.  Using 
dopants  of  different  density  or  atomic  number  than  diose  given  above  has  also  been  considered,  but 
detailed  calculations  have  shown  such  approaches  to  be  less  dian  optimum.  For  very  low  density 
foams  in  which  doping  is  not  practical,  die  use  of  a  foam  contiuning  oxygen  is  recommended  as 
our  diow  the  oxygen  lines  may  potentially  contain  diagnostic  information  which 

is  difficult  to  obtain  from  carbon  alone. 

We  have  also  carried  out  a  full  radiation-hydrottynamic  emulation  of  an  experiment  at 
Sandia  Natitmal  Laboratories  in  which  a  Z  pinch  x-n^  source  heated  a  1  /im  diick  **macrodot’’  of 
MgFj.  The  MgF]  gradually  heats  and  esqiands  and  near  die  peak  of  the  radiation  pulse  reaches 
tonperatures  of 70-80  cV  and  electron  densities  of  (1-4)  x  10*®  cm"®.  These  values  are  consistent 
with  previous  itgrimiitits  at  Sandia  using  line  ratios,  indicating  diat  a  good  understanding  of  die 
experiment  has  been  achieved. 


E-1 


TEMPERATURE  DIAGNOSIS  OF  lON-BEAM-DRIVEN  FOAM  TARGETS 
USING  K  SHELL  FEATURES  FROM  TRACE  ELEMENTS 


I.  Introduction 

In  evaluating  the  performance  of  pulsed-power-driven  targets,  the  need  to  determine  their 
temperature  frequently  arises.  This  is  particularly  true  for  foam  converters  used  in  inerrial 
confinement  fusion  (ICF)  driven  by  light  ion  beams.  At  Sandia  National  Laboratories,  the  Particle 
Beam  Fusion  Accelerator  n  (PBFA 11)  deposits  ion  energy  into  a  low-density  hydrocarbon  foam 
target  for  the  purpose  of  converting  this  energy  to  copious  x  radiation  to  drive  a  fuel  crq>sule.  This 
indirect  drive  method,  aimed  ultimately  at  achieving  breakeven  and  beyond,  requires  nprimnm  and 
well  understood  performance  from  the  foam  plasma  radiation  converts. 

A  key  parameter  in  determining  the  radiation  temperature  and  spectrum  is  the  electron 
temperature  of  the  ion-beam-driven  foam  plasma.  In  addition  to  being  one  of  the  benchmarks  of 
converter  performance,  the  foam  plasma  temperature  is  also  an  important  indicator  of  how  well  the 
ion  beam  has  coupled  to  the  target.  Many  of  diese  issues  are  discussed  in  greater  detail  elsewhere^. 
As  indicated  in  Ref.  1 ,  many  diagnostic  instruments  have  been  successfully  dq>loyed  in  the  difficult 
environment  of  the  PBFA  n  pulsed  ion  accelerator,  among  them  several  spectrogrrq^  and  x-ray 
diodes.  However,  some  useful  diagnostic  techniques,  such  as  laser-driven  target  haclrlighting,  are 
not  currently  feasible  in  PBFA  n  experiments  due  to  cost  and  ^ace  limitations.  The  purpose  of  this 
report  is  to  assess  in  detail  a  flexible  yet  simple  and  accurate  method  that  m  akes  use  of  conventional 
x-ray  q>ectroscopy  and  is  practical  for  PBFA  n  targets.  In  this  iq>proach,  the  hydrocarbon  foam 
converters  are  doped  with  trace  amounts  of  elemoits  (O,  F,  Na)  \tdiich  absorb  or  radiate  different 
K  shell  features  from  different  ionization  stages  as  the  temperature  changes.  We  emiriiasize  that 
the  presently  described  techniques  i^ly  to  the  target  after  the  passage  of  the  energizing  ion  beam, 
therefore  the  spectral  responses  due  to  direct  excitation  of  atomic  levels  by  the  ion  beam  ate  not 
modeled. 

In  the  succeeding  sections  this  technique  is  evaluated  in  detail  for  a  variety  of  expected  target 
conditions.  Sec.  n  describes  the  NRL  model  for  coupled  atomic  physics  and  radiation  in  aplasma, 
and  compares  our  calculated  opacities  to  both  the  XSN  and  Los  Alamos  values  for  parameters 
relevant  to  expected  PBFA  n  conditions.  In  Sec.  m  we  discuss  the  basic  radiation  j^ysics  which 
underlies  the  formation  of  the  principal  x-ray  features,  whether  emission  or  absorptitm.  Sec.  IV 
presents  results  for  homogeneous  doped  foams  with  both  temperature  and  density  varying  or  with 
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fixed  density  and  variable  temperature.  Sec.  V  discusses  inhomogeneous  foam  plasmas,  and  in 
Sec.  VI,  a  full  radiation-hydrodynamic  simulation  of  the  MgFj  “macrodot**  e3q)etiment  carried  out 
on  the  Sandia  pulsed  power  accelerator,  Saturn,  is  presented. 

n.  Model  and  Opacity  Comparison 

The  results  presented  below  were  obtained  with  a  plasma  model  which  has  been  in  use  for 
some  time  at  the  Naval  Research  Laboratory  and  is  extensively  documented  in  tiie  literature. 
This  model  emplo3rs  detailed  configuratitm  structures  of  various  atomic  ions  and  tiieir  associated 
coupling  rates  and  solves  for  the  self-omsistent  le>^  populations  and  radiation  field  (botii  internal 
and  emitted)  in  a  ^)ecified  plasma.  Radiative  transfer  is  calculated  nring  either  a  cdl-to-cell 
coupled  escape  probability  technique  or  full  multiftequency  treatment.  Due  to  die  significant  line- 
ctmtinuum  interactions  in  tile  ion  beam  taiget  plasmas  analyzed  in  die  present  wodt,  multifirequency 
tranqwft  was  emplt^ed.  Various  assumptions  can  be  invdced  in  exercising  die  modd,  namdy 
Local  Tbermocfynamic  Equilibrium  (LTE),  Collisional  Radiative  Equilibrium  (C!RE),  or  die  full 
time-dependent  rate  equations  can  be  integrated.  The  latter,  obviously  mcneeaqiensiveqjproach  is 
usually  reserved  for  cases  where  die  atmnic  {d^cs  and  radiation  are  couided  to  a  hydrodynamic 
simulatitm  (see  Ref.  2,  for  example).  The  foil  multifrequencytran^xm  option  has  been  employed 
in  previous  instances  where  ICF  targets  (laser-driven  cqisules)  were  investigated.  The  NRL 
model  obtained  excellent  agreement  with  experimentally  observed  x-n^  qiectra  for  implosims  of 
neon-filled  microballoons  at  the  Univeisity  of  Roctester  (Ref.  3)  and  DT-argon  filled  qdierical 
targets  at  Los  Alamos  Naticmal  Laboratory  (Ref.  4). 

In  the  targets  considered  in  this  work  several  elements  ate  present.  In  addition  to  dw  carbon 
and  hydrogen  of  die  foam,  dopants  (O,  F,  and  Na)  are  used  as  tracer  temperature  diagnostics.  A 
separate  data  table  is  consulted  for  each  element;  energies  of  levels  and  radiative  rates  are  obtained 
from  various  literature  sources  or  Cowan’s  atmnic  structure  code.  Collisioiial  and  {diotoitniization 
rates  are  fonn  distorted  wave  or  semidassical  calculations,  or  previoudy  published  sources,  as 
appropriate.  We  aim  to  include  all  important  atmnic  processes  coupling  the  modeled  levels,  such  as 
dectron  collisional  mcdtatitm  and  de-exdtatitm,  radiative,  3-body  and  didectrmuc  recombinatirm, 
coUisitHial  and  jdiotoitniizatimi,  bound-bound  radiative  decay  and  photoexdtation.  In  solving  the 
radiative  transfer  equatirm,  the  qpadty  at  each  of  qiproximately  2600  freque-ides  is  obtained  by 
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adding  the  absorption  coefficients  from  bound-bound,  bound-free,  and  free-free  radiative  processes 
for  each  of  the  elements  in  the  plasma.  This  generates  the  internal  radiation  field  for  calculation 
of  photoexdtation  and  photoionization  rates  as  well  as  the  emitted  spectrum  for  comparison  with 
experimental  data.  A  Voigt  profile  is  assumed  for  each  of  the  spectral  lines,  with  the  broadening 
parameter  obtained  from  the  radiative  and  coUisional  lifetimes  of  the  upper  and  lower  levels  of  the 
transition  in  question.  Further  details  regarding  the  rates  may  be  found  in  Ref.  5. 

For  the  dopants,  O,  F,  and  Na,  as  well  as  Mg  which  is  needed  for  a  separate  simulation,  the 
atomic  models  used  are  analogous  for  isoelectronic  sequences.  Exdted  state  structure  is  induded 
for  the  lithiumlike,  heliumlike,  and  hydrogenic  ions  through  n=5.  The  hydrogenic  levels  are 
differentiated  by  prindpal  quantum  number  only.  For  the  heliumlike  stage,  the  structure  is  broken 
into  individual  singlet  and  triplet  sublevels  throu^  n=3,  and  treated  as  single  Rydberg  levels  by 
prindpal  quantum  number  only  for  higher-lying  states.  In  die  lithiiimlike  sequence,  the  levels  are 
brdcen  out  by  bodi  prindpal  and  orbital  quantum  number  through  n=3,  and  for  na4  and  5,  tmly 
by  prindpal  quantum  number.  For  carbon  a  somewhat  simpler  model  is  sometimes  used  which 
includes  states  only  through  n=2  for  the  K  shell.  This  reduces  cost  and  storage,  and  does  not 
sacrifice  essential  accuracy  because  the  electrcm  densities  considered  are  generally  suffident  to 
force  LTE  pt^iulations  for  carbon  but  not  for  the  higher  Z  dopants,  as  discussed  below  in  the  next 
section.  The  accurades  of  die  qpadties  given  by  this  carbon  model  are  demonstrated  bdow  in  a 
comparison  widi  other  models.  Note  that  the  hydrogen  is  nearly  always  fully  ioniz^  and  does 
litde  other  than  contribute  electrons  and  add  some  mass.  It  is,  however,  modeled  with  a  small 
atomic  data  table  containing  the  neutral  ground  and  bare  nudear  stages. 

Our  model  generates  opadties  comparable  to  those  obtainable  frtHn  two  widely  employed 
data  bases,  the  XSN  model*,  developed  at  Lawroice  Livermore  National  Laboratory,  and  the  Los 
Alamos  opadties^ .  We  have  chosen  conditions  for  comparison  which  are  typical  of  those  used  in 
our  foam  model  calculations.  The  mass  density  is  20  mg  cm~*  of  CH,  corresponding  to  an  attnnic 
density  of  9.3  x  10^*  cm~*.  The  electron  temperature  is  lOOeV.  Fig.  1  presents  a  direct  comparison 
of  the  detailed  spectral  opadties  predicted  by  the  three  models.  Note  the  genoral  agreement  across 
the  spectrum  extending  from  10  eV  to  10  keV.  Minor  differences  are  largely  due  to  the  presence 
or  absence  of  specific  lines  and  their  assumed  profiles.  Note  that  the  opadty  is  givoi  in  cm^  g~^ , 
easily  convertible  to  absorption  coeffident  in  cm~^  by  multiplying  by  the  density,  in  diis  case  0.02 
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Efwrsy  (•V) 


no.  1.  OpadQfiiicm^  g"*  ifploctedagaiiisteneigyfiom  lOeVtolOlceVforCHatandectron 
tempentture  of  100  cV  md  num  denai^  of  20  mg  cm"*.  Hieiiesoltsof  Ane^fEefeatmodds, 
NRL,  XSN,  and  LANL,  are  diown  for  oompaiuon. 
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g  cm~®.  At  energies  below  100  eV,  free-free  absorption  dominates  the  opacity,  whereas  above 
this  value,  bound-free  transitions  are  the  principal  contributor.  In  oxygen-free  foams,  carbon  is 
die  principal  source  of  opacity,  and  this  comparison  demonstrates  that  our  model  can  accurately 
characterize  this  critical  parameter  in  the  regime  of  interest. 

m.  Relevant  Radiation  Physics 

The  purpo%  of  this  section  is  to  describe  the  jdiysical  principles  which  explain  die  qiecific 
results  presmted  in  the  succeeding  sections.  Those  results  suggest  observable  qiectroscopic 
features  as  diagnostic  tools  for  die  amditions  of  die  target  plasma.  Therefore  we  ask  die  basic 
question:  what  governs  the  formation  of  x-n^  spectra  in  the  types  of  target  plasmas  of  interest? 

For  a  plaana  to  radiate  as  a  blackbody  two  conditions  must  be  met.  First,  the  plasma  must  be 
optically  thick  in  the  spectral  region  being  considered,  and,  second,  the  level  populations  which 
do  die  radiating  must  be  distributed  according  to  die  Saha-Boltzmaim  equation.  The  first  of 
these  conditions  is  straightforward  and  easily  understood.  One  simply  multiplies  die  absorption 
coefficient  by  the  thickness  of  die  plasma.  We  consider  foam  plasmas  0.5  cm  diick  for  most  of  die 
cases  run  in  diis  study.  The  optical  dqidifortheconditimis  of  Hg.  1,  for  example,  may  be  obtained 
by  multiplying  the  opacity  in  cm^  g~^  by  the  density  (0.02  g  cm~’)  and  die  0.5  cm  thickness, 
for  any  jdioton  oieigy  of  interest.  When  dus  is  done  for  Fig.  1,  it  is  found  that  die  optical  dqidi 
exceeds  unity  excqit  for  a  small  region  of  energies  near  4(X)  eV,  and  even  there,  die  optical  depdi 
is  a  few  tenths  or  so.  Therefore,  a  CH  foam  of  demsity  20  mg  cm~’  has  the  potential  to  radiate 
a  near-Planddan  spectrum  over  a  large  regicm  of  the  soft  x-ray  spectrum,  if  its  populaticms  are  in 
LTE,  which  we  now  consider. 

For  a  Saha-Boltzmarm  (LTE)  prqiulation  distribution  to  prevail,  die  partide  distribution  must 
be  Maxwellian,  which  we  assume  throughout  tins  work,  and  the  collision  rates  produced  by  those 
particles  must  well  exceed  the  corre^ionding  radiative  decay  rates.  We  confine  our  discussion  to 
the  K  shell  ionization  stages,  for  which  die  spontaneous  decay  rates  vary  as  (Ref.  8),  where 
Z  is  the  atomic  number.  The  collisional  rates  vary  as  NeZ~’  (Ref.  9),  therefore  to  maintain  a 
fixed  ratio  of  collisional  to  radiative  rates  required  for  LTE  means  that  the  electron  density  must 
scale  as  7? .  Qearly,  it  is  much  easier  to  bring  the  K  shell  into  LTE  for  low  atomic  number 
elements.  The  detailed  investigation  of  Ref.  10  found  that  for  optically  thin  A1  (Z  s  13),  an 
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electron  density  of  about  10^®  is  needed  to  force  the  K  sheU  stages  into  LTE.  The  electron  density 
required  for  LTE,  assuming  7}  scaling,  is  plotted  in  Fig.  2  for  atomic  numbers  6  <  Z  <  13.  Note 
that  LTE  for  K  shell  carbon  requires  an  electron  dmsity  of  ~  5  x  10’’  cm~’,  about  an  order  of 
magnitude  higher  than  that  obtained  at  a  foam  density  of  20  mg  cm~’.  However,  virtually  all  of 
the  principal  resonance  transitions  are  optically  thick  at  diis  density  for  a  half  centimeter  foam. 
The  resulting  self-absorption  in  the  lines  greatly  reduces  the  effective  decay  rate,  thus  the  electron 
(tensity  iweded  for  LTE  populations  is  effectively  lowered  to  values  close  to  tfiose  attained.  It  is 
therefore  reasonable  to  expect  cartxm  to  be  faiiiy  dose  to  LTE.  However,  Na,  one  of  die  dopants 
considered,  requires  (tensity  ~  (11/6)’  ^  70  times  that  of  C  and  dierefore  should  not  be  expected 
to  be  nearly  as  close  to  LTE. 

To  illustrate  this  point  and  also  show  why  absorption  lines  are  prevalent  in  much  of  the 
results  that  follow,  consider  the  excitadon  of  the  principal  heliumlike  resonance  line  of  Na  at  its 
energy  E  of  1127  eV.  This  is  a  very  strong  transititm  with  an  oscillator  strength,  /,  of  0.75.  For 
idiotoexdtation,  assume  that  the  sodium  ion  is  in  the  interior  of  the  foam  and  sees  a  100  eV 
Planddan  spectrum  over  all  Ak  steradians.  The  photoexdtadon  rate  P  is  given  by 

/.(„.-■)  =  I  (1) 


where  <r„  is  the  line  cross  sectitm  in  cm’,  whose  fretpiency  integral  is  ^/.  Using  die  excellent 
^^roximation  that  the  Planck  function  is  constant  over  the  narrow  spectral  tini»  profile  leads  to 


P  = 


4.34  X  10’£?’(eF)  , 

eE/kT  _i  i' 


(2) 


Eq.  (2)  yields  a  line  jdiotoexdtation  rate  of  5.27  X 10*  sec~^  for  the  s(xlium  He  a  line  immersed  in 
a  100eVblackbo(ty.  At  the  surface  of  the  foam  plasma,  the  ambient  radiation  fidd  is  reduced  by  at 
least  a  factor  of  2  due  to  geometry  al(me.  Greater  reductions  are  possible  due  to  line  self-absorpti(m 
(see  Ref.  11).  By  contrast,  the  collisional  exdtation  rate  of  the  upper  level  of  this  line  at  100 
eV  and  electron  density  of  5  x  10’^  cm“’  is  just  6  x  10*  sec“^  Accordingly,  it  is  expected 
that  the  upper  level  populations  of  the  dopants  may  be  dominated  by  photoexdtation  and  dius  be 
reduced  near  dw  boundaries  of  the  foam.  Since  many  of  their  K  shell  lines  are  optically  thirlc^  this 
has  the  consequence  that  they  may  be  seen  in  absorption  against  a  background  of  near-Planddan 
continuum  radiation.  This  is  predicted  to  be  the  case  (]uite  often,  as  seen  in  the  following  sectitm. 
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IV.  Homogeneous  Doped  Hydrocarbon  Foam  Plasmas 

The  following  subsections  report  the  results  of  detailed  calculations  in  which  trace  amounts  of 
elements  of  similar  atomic  number  are  assumed  doped  into  the  foam  target.  A  similar  technique  has 
been  used  with  considerable  success  at  the  University  of  Rochester^^  to  diagnose  the  temperature 
of  laser-imploded  spherical  targets.  The  advantages  of  this  method  are  discussed  in  Ref.  12.  An 
important  point  is  that  the  functional  forms  of  the  tonperature  and  atomic  number  dq>endences 
of  the  various  rate  coefficients  are  well  known.  Even  thougfi  the  absolute  values  of  the  rates 
may  contain  significant  errors  and  uncertainties,  in  taking  their  ratios  from  two  elements  of 
nearby  atomic  number,  those  errors  generally  cancel  to  first  order.  This  is  especially  true  for  the 
well-known  atomic  physics  of  the  K  shell,  to  vdiich  we  confine  our  analysis.  In  the  temperature 
regime  of  interest  for  PBFA  n  experiments,  elements  with  atomic  number  above  sodium  (Z=l  1) 
are  only  weakly  excited  in  die  K  shell.  The  most  promising  diagnostic  tracers  were  found  to  be  F 
and  Na  (which  presumably  can  be  easily  doped  as  NaF)  as  well  as  O,  which  can  also  be  present  as 
integral  to  the  foam  target  composition.  Necm  would  probably  also  be  excellent  but  its  status  as  a 
noble  gas  in  all  likelihood  renders  its  use  unfeasible  and  it  was  not  considered  in  the  present  work. 

A.  Temperature  Increase  with  Density 

The  initial  studies  crmsider  foams  in  which  die  density  and  temperature  are  assumed  to 
increase  concunendy.  Four  specific  CH  foams  were  analyzed  in  detail:  assumed  temperatures  and 
densities  of  50  eV  and  10  mg  cm"’,  1(X)  eV  and  20  mg  cm"’,  150  eV  and  30  mg  cm"’,  and  a 
200  eV  foam  with  an  assumed  density  of  40  mg  cm"’.  Clalculations  were  run  in  which  uniform 
doping  with  NaF  at  various  densities  was  assumed,  and  also  in  which  a  solid-density  layer  of  NaF 
was  placed  on  one  side  of  the  foam. 

The  most  promising  configuration  was  found  to  be  uniform  doping  at  a  concentration  of  6  x 
10^’  atoms  cm"’  of  both  Na  and  F.  An  extensive  series  of  computations  at  lower  densities  showed 
that  no  additional  diagnostic  insights  were  obtained,  but  the  characteristic  features  were  weakened 
and  therefore  more  difficult  to  detect.  At  6  x  10^®  cm"’ ,  the  mass  density  of  NaF  is  0.4  mg  an"’ , 
unlikely  to  affea  the  hydrodynamics  of  the  foams  which  are  well  over  an  order  of  magnitude  more 
dense.  However,  at  an  order  of  magnitude  higher  dopant  concentrations,  the  NaF  mass  density 
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would  be  an  appreciable  fraction  of  the  foam  density  and  thus  perturb  its  behavior. 

In  calculations  in  which  a  contiguous  solid  density  layer  of  NaF  is  present,  the  electron 
density  of  about  3  x  10^’  cm~^  greatly  enhances  the  3-body  recombination  rate,  making  it  much 
harder  to  ionize  these  elements  to  the  K  shell.  These  two  elements  are  also  much  closer  to  LTE 
than  at  dopant  densities.  Therefore,  there  is  a  nearly  LTE  solid  layer  radiating  against  a  nearly 
Planddan  continuum  produced  by  the  carbon  in  the  foam.  Not  surprisingly,  the  features  which  are 
present  are  very  weak  and  of  little  diagnostic  value.  The  specific  results  of  the  most  promising 
configuration,  namely,  the  uniform  NaF  doping  at  6  x  10^‘  cm~^ ,  are  now  considered. 

Hgures  3-fi  |msent  tiie  region  of  tiie  plectrum  from  0.7  to  1 .7  keV,  as  calculated  to  be  emitted 
by  an  0.5  cm  thick  foam  at  temperatures  of  50, 100, 150,  and  200  eV,  at  CH  densities  of  10, 20, 
30,  and  40  mg  cm“®,  respectively.  Uniform  dt^nng  of  NaF  at  a  density  of  6  x  10^®  an~®  is 
assumed  in  each  case.  One  quite  attractive  feature  of  using  NaF  is  that  their  K  shell  lines  are 
cleanly  separated,  all  of  the  F  lines  lying  below  1100  eV  and  the  Na  lines  being  above  tiiis  energy. 
Note  from  Figs.  3-6  that  most  of  the  lines  appear  in  absoiption,  for  reasons  discussed  above  in 
Sec.  m.  The  results  show  that  the  various  temperatures  can  be  distinguished  experimentally  by 
mostly  obvious  qualitative  features  of  the  spectra.  In  each  figure  the  dashed  line  is  the  blackbody 
spectrum  for  the  calculated  temperature  shown  for  comparison.  The  distinguishing  characteristics: 

•  50  eV:  Only  the  He-like  lines  of  F  are  evident,  all  in  absorption. 

•  100  eV:  Both  the  He-and  H-like  lines  of  F  are  seen  in  absorption;  for  Na,  only  the  He-like 
lines  are  seen,  all  in  absorption. 

•  150  eV:  Both  He-and  H-Uke  lines  of  Na  are  now  seen  in  absorption:  both  species  of  F  are 
still  visible.  A  hint  of  emission  sppears  in  the  wings  of  some  of  the  higher  lines  in  Na. 

•  200  eV:  The  Lyman  lines  of  F  arising  from  n  >  2  and  all  lines  of  Na  for  n  >  2  are  now  seen 
in  emission. 

These  calculations  establish  clear  benchmarks  for  determining  the  foam  temperature  with  ~ 
50  eV  resolution  from  50  to  200  eV,  when  the  density  increases  along  with  the  temperature.  It  is 
of  interest  to  also  determine  how  well  one  can  rejolve  the  temperature  spectroscopically  with  this 
technique,  and  also  if  increasing  the  temperature  at  a  fixed  density  leads  to  significant  quantitative 


9 
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FIO.  3.  Emitted  spectrum  ss  calculsted  to  srise  from  s  CH  fosm  plssms  of  density  10  mg[  cm~* 
and  temperature  50  eV.  The  thickness  of  the  foam  is  0.5  cm  and  it  is  assumed  to  be  doped  with  NaF 
at  a  density  of  6  x  10^*  cm”* .  The  dashed  line  is  a  50  eV  blackbody  spectrum  for  comparison. 
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FIG.  4.  As  in  Fig.  3,  except  that  the  plasma 
respectively. 
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or  qualitative  changes. 


B.  Fixed  Density  Calculations 

Tte  spectra  discussed  in  this  section  were  all  calculated  assuming  a  fixed  CH  foam  density 
of  10  mg  cm~^.  All  dopants  are  spatially  uniform  at  an  atomic  density  of  6  x  10^’  cm~’.  The 
temperatures  considered  are:  SO,  60,  70,  80, 90, 100, 120, 140, 150,  and  200  eV.  The  req)ective 
spectra  appear  in  Figs.  7-14. 

Initial  calculations  were  petfonned  with  NaF  as  die  only  dopant,  as  was  done  in  die  woik 
discussed  in  the  previous  section.  Clear  and  useful  diagnostic  trends  were  seen  widi  promising 
experimental  boichmaiks  for  each  increment  in  temperature,  excqtt  for  die  60-70  eV  interval, 
where  only  mild  sensitivity  was  observed  in  the  Na  and  F  spectra.  Therefore,  an  oxygen  K  shell 
table  was  set  up,  and  O  was  assumed  as  adtqiant  (also  at  an  atomic  density  of  6  x  10^*  cm~’)  for 
therunsfirQmS0to70eV.  For  these  three  q)ectra,idioton  energies  from  500  to  15(X)  eV  are  shown 
so  diat  the  oxygen  K  shell  lines  may  be  included  in  the  plots,  whereas  the  odier  runs  only  use  Na 
and  F  and  the  spectra  are  displayed  from  700  to  1700  eV.  Summarizing  the  important  features: 

•  50eV:  Only  the  He-like  lines  of  both  F  and  O  are  present,  all  in  absorption. 

•  60  eV:  The  He-like  lines  of  Na  now  qipear. 

•  70  eV:  The  Ly  a  line  of  O  qipears. 

•  80  eV:  The  Ly  a  line  of  F  iqipears. 

•  90  eV:  The  higher  Lyman  series  lines  of  F  are  now  detectable. 

•  100  eV:  As  the  temperature  increases,  the  continuum  begins  to  fall  away  fin»n  the  Planck 
limit  in  the  keV  region.  This  causes  emission  lines  to  form  for  the  higher  series  He-like  Na  lines. 
The  density  is  half  of  the  20  mg  cm~’  assumed  in  die  previous  calculation  at  100  eV.  The  emission 
lines  have  self-reversed  cores. 

•  120  eV:  The  F  lines  arising  frxim  n  >  3  are  now  in  emission. 

•  140  eV:  F  Ly  a  is  now  in  emission  but  with  a  self-reversed  core. 

•  150  eV:  The  Ly  a  line  of  Na  qipears  with  a  deep  self-reversed  core. 


Intensity  (ergs/sec-cm^-keV) 


Intensity  (ergs/sec-cm^-keV) 
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Intensity  (ergs/sec-cm^-keV) 


Intensity  (crgs/sec-cm^-keV) 


Inlensily  (ergs/sec-cm^-keV) 


nC.  14.  The  predicted  K  shell  spectrum  ofdie  sodium  dopant  changes  markedly  as  the  temperature 
increases  from  100  to  200  eV.  Foam  density  is  fixed  at  10  mg  cm~’ . 


•  200  eV:  The  foam  density  of  10  mg  cm~^  is  only  one-fourth  of  that  used  in  the  previous 
series  of  calculations.  As  a  result,  the  continuum  is  so  far  below  the  Planck  limit  that  all  the  lin«»!8 
are  now  in  emission. 

The  principal  conclusions  of  our  study  at  a  fixed  density  of  10  mg  cm~^  are  :  i)  temperature 
distinctions  as  fine  as  10  eV  may  be  made  in  the  50-100  eV  regime  i»oviding  that  oxygen  is 
included  as  a  dopant  along  with  Na  and  F,  ii)  the  fiu(»ine  spectrum  can  differentiate  die  temperature 
in  20  eV  increments  ffom  80  eV  to  140  eV,  and,  iii)  due  to  the  lower  carbon  density,  the  continuum 
falls  further  firom  tiie  Planck  limit,  especially  at  the  higher  temperatures  cmsidered.  The  result  is 
that  most  of  the  spectral  lines  are  seen  in  emission  at  150  eV,  aixl  all  at  200  eV. 

V.  Inhomogeneous  Foam  Plasmas 

In  tile  initial  ion-beam  foam  target  expetimoits  tiie  foams  expected  to  be  used  are  of  somewhat 
lower  density  than  those  considered  above  (i.e.,  about  3  mg  cm~’  instead  of  10  mg  cm~’  or 
higher).  Doping  of  these  foams  with  tracer  elements  is  difficult  at  such  low  densities  although 
it  is  expected  to  be  feasible  eventually,  and  is  certainly  possible  at  hi^ier  target  densities.  To 
provide  initial  diagnostic  piedictimis  for  tiiese  experiments  we  have  performed  calculations  for  the 
emitted  spectra  of  foams  of  CH  as  well  as  CuHisOio.  Ihe  temperature-density  profile  is  a  more 
realistic  non-uniform  distribution  from  a  radiaticm-hydrodynamics  calculation,  courtesy  of  M.  K. 
Matzen  and  R.  Dukart  of  Sandia  National  Laboratories.  The  same  profile  is  used  for  both  foam 
compositions  and  is  illustrated  in  Hgs.  15  and  16. 

The  q)ectra  presented  in  Figs.  17-19  following  the  temperature  and  density  plots  are 
calculated  as  viewed  from  the  low-temperatuie  side  of  the  foam  (displacement  s  0).  Fig.  17 
shows  the  entire  spectrum  of  die  CH  foam.  A  somewhat  more  detailed  carbcxi  table  has  been 
used  than  that  employed  for  previous  calculations,  allowing  the  higher  level  K  shell  lines  to  be 
included.  Note  from  Fig.  17  that  the  lower  plasma  density  has  resulted  in  the  spectrum  falling 
significantly  below  the  Planck  limit,  although  its  overall  shape  could  still  be  credibly  interpreted 
as  Planddan.  The  comparison  blackbody  curve  is  for  an  85  eV  radiatirm  temperature  since  tiiat 
is  the  highest  temperature  attained  in  the  foam  plasma.  Fig.  18  shows  a  closeup  of  the  carbon 
K  shell  lines,  with  the  Ly  a  -  He  7  blend  the  most  prominent  feature,  and  He  a  nearly  of  equal 
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displacement  (cm) 

FIG.  15.  Calculated  temperature  profile  of  ion  beam  foam  target  converter,  courtesy  of  M.  K. 
Matzen  and  R.  Dukart  of  Sandia  National  Laboratories. 
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PIO.  18.  Portion  of  the  same  spectnun  of  Fig.  17,  focusing  on  the  K  shell  lines  of  carbon. 
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intensity.  Fig.  19  illustrates  the  spectrum  from  an  oxygen-bearing  foam  (C12H18O10)  of  the 
identical  temperature-density  profile  as  the  CH  case.  All  of  the  principal  K  shell  lines  of  oxygoi 
appeal  in  absorption  against  the  strong  hydrogoiic  continuum  emission  of  carbon.  As  with  sodium 
and  fluorine,  the  carbon  and  oxygen  K  lines  are  cleanly  separated,  all  of  the  carbon  lines  lying 
below  500  eV,  and  the  oxygen  lines  all  above  SSO  eV.  The  strong  presence  of  O  Ly  a  indicates  that 
temperatures  of  at  least  70  eV  have  been  attained  in  the  target.  This  temperature  benchmaik  using 
oxygen  has  also  beat  noted  above  for  previous  calculations.  Therefore,  even  if  dopants  caiuiot  be 
used  in  a  particular  experiment,  use  of  oxygen-bearing  foams  is  recommended  as  a  way  to  obtain 
additional  diagnostic  infonnation. 

VI.  Simulation  of  Magnesium  Fluoride  Macrodot  Experiment 

In  a  recent  experiment  at  Sandia,  a  Ifon  tfaidc  target  of  MgFj  (referred  to  as  a  **macrodot**) 
was  exposed  to  the  x-rays  from  a  Z  jnndi  radiation  source.  The  reqxmse  of  the  MgF}  target 
to  the  x-rays  was  observed  with  x-ray  spectrometers^’,  which  were  looking  primarily  at  the  Mg 
K  shell  lines.  In  this  spectral  region,  the  He  a  liiM  was  the  dominant  feature,  and  fire  He-like 
intercombination  line  and  various  satellite  lines  were  also  detected.  An  initial  analysis^’  employed 
standard  line  ratios  involving  the  restmance,  inteicombmation,  and  satellite  lines.  This  mediod 
suggested  that  temperatures  of  50-70  eV  had  been  attained,  at  electron  densities  of  10’°  -  10’^ 
cm“’. 

Using  the  multifrequency,  one-dimensional  version  of  the  NRL  radiation  hydrodynanrics 
model,  we  have  simulated  this  experiment  widi  the  atomic  data  tables  for  Mg  and  F  described 
above  in  Sec.  n.  The  Z  pinch  source  spectrum  is  modeled  as  a  two-componem  radiation  field.  Both 
components  have  a  Gaussian  time  profile,  with  a  full-width-at-half-maximum  of  15  ns.  The  low 
energy  ("soft”)  componoit  is  taken  as  an  undiluted  48  eV  blackbody  at  the  peak  of  the  Gaussian. 
The  temperature  of  the  assumed  Planckian  varies  such  that  T*  is  Gaussian.  The  higher  energy 
(“hard”)  component  is  assumed  to  consist  of  photons  lying  entirely  between  1.6  and  2.1  keV, 
whose  intensity  per  unit  energy  interval  is  constant  in  that  range.  The  total  hard  power  is  120  GW 
cm~’  at  the  peak  of  the  radiation  pulse. 

The  calculation  reveals  that  the  l/^m  of  MgF2  heats  and  expands,  much  like  exploding  foils 
used  in  x-ray  laser  experiments^^  or  the  plastic  window  used  to  retain  the  neon  in  the  Sandia 
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sodium-neon  x-ray  laser  targets^®.  Fig.  20  shows  that  a  gentle  temperature  gradient  is  established 
by  the  time  of  the  peak  of  the  external  x-ray  pulse.  The  calculated  temperature  of  70-80  eV  is  in 
satisfactory  agreement  with  the  initial  estimates^®  of  50-70  eV.  Furthermore,  the  electron  density 
as  shown  in  Fig.  21  varies  from  (1-4)  x  10^“  cm“®,  also  in  good  agreement  with  the  Ref.  13 
analysis.  Fig.  22  shows  the  calculated  spectrum  from  0.7  to  2.0  keV,  containing  all  the  K  shell 
lines  of  both  F  and  Mg.  This  spectrum  reveals  that  Mg  He  a  is  predicted  to  be  the  highest  intensity 
line,  also  in  agreemmt  with  experiment.  This  very  respectable  convergence  on  conditions  using 
indep«ident  methods  demonstrates  a  good  understanding  of  the  physics  of  the  experiment  as  well 
as  a  useful  benchmark  for  the  model. 

Finally,  we  note  that  two  fundamental  processes  are  responsible  for  stripping  the  Mg  plasma. 
One  is  heating  of  the  plasma  by  the  radiation  source,  which  makes  the  electrons  hot  enough  to 
further  coUisionally  ionize  the  Mg.  The  other  is  direct  photoionization  of  the  Mg  ions  die 
radiation  field.  The  effect  of  direct  {diotoionization  is  shown  in  the  calculations  presented  in  Figs. 
23  and  24,  where  the  He-like  fraction  is  shown  as  a  function  of  temperature  and  density.  The 
contours  for  a  He-like  fraction  exceeding  0.9  are  considerably  broader  in  Fig.  24,  where  the  effect 
of  the  pump  field  is  included,  than  in  Fig.  23,  whm  it  is  not.  This  effect  is  most  pronounced  at 
the  lower  densities  as  would  be  expected,  since  colMonal  ionization  is  weaker. 
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FIG.  20.  Temperature  of  the  MgFa  *^acrodot”  plasma  vs.  position  at  the  peak  of  the  external 
heating  pulse.  Initial  position  of  the  target  is  diq>lacement «  0. 
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ke  Fractional  Population 


Electron  Temperature  (eV) 

FIG,  23.  He-like  fraction  for  Mg  is  plotted  as  a  function  of  temperature  and  density  with  n( 
external  radiation  field  pumping  the  plasma  in  the  calculation.  Collisional-radiative  equilibrium  i 
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Electron  Temperature  (eV) 

FIG.  24.  As  in  Fig.  23,  except  that  photoionization  by  the  external  pump  radiation  source  discussed 
in  Sec.  VI  is  included  in  calculating  the  ionization  balance. 
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